Small-scale clumps in the galactic halo and dark matter annihilation 
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Production of small-scale DM clumps is studied in the standard cosmological scenario with an 
inflation-produced primeval fluctuation spectrum. Special attention is given to the three following 
problems: (i) The mass spectrum of small-scale clumps with M < 10 3 Mq is calculated with tidal 
destruction of the clumps taken into account within the hierarchical model of clump structure. Only 
0.1 - 0.5% of small clumps survive the stage of tidal destruction in each logarithmic mass interval 
AlnAf ~ 1. (ii) The mass distribution of clumps has a cutoff at M m i n due to diffusion of DM 
particles out of a fluctuation and free streaming at later stage. M m in is a model dependent quantity. 
In the case the neutralino, considered as a pure bino, is a DM particle, M min ~ lO~ 8 Af . (iii) The 
evolution of density profile in a DM clump does not result in the singularity because of formation of 
the core under influence of tidal interaction. The radius of the core is R c ~ 0.1R, where R is radius 
of the clump. The applications for annihilation of DM particles in the Galactic halo are studied. The 
number density of clumps as a function of their mass, radius and distance to the Galactic center 
is presented. The enhancement of annihilation signal due to dumpiness, valid for arbitrary DM 
particles, is calculated. In spite of small survival probability, the global annihilation signal in most 
cases is dominated by clumps. For observationally preferable value of index of primeval fluctuation 
spectrum n p « 1, the enhancement of annihilation signal is described by factor 2-5 for different 
density profiles in a clump. 

PACS numbers: 12.60.Jv, 95.35.+d, 98.35. Gi 



I. INTRODUCTION 

Both analytic calculations PJ Q and numerical simu- 
lations 0. [J Q predict existence of dark matter clumps 
in the Galactic halo. The density profile in these clumps 
according to analytic calculations [H 0, H E| an d nu- 
merical simulations 0, ^} is p(r) oc . An average 
density of the dark matter (DM) in a galactic halo itself 
also exhibits a similar density profile (relative to a galac- 
tic center) in the both approaches. The DM profiles in 
clusters of galaxies is discussed in ^lj an d in references 
therein. In the analytic approach of Gurevich and Zybin 
(see review and references therein) the density pro- 
files are predicted to be universal, with (3 « 1.7 — 1.9 for 
clumps, galaxies and two-point correlation functions of 
galaxies. In the numerical simulations the density pro- 
files can be evaluated only for the relatively large scales 
due to the limited mass resolution. The value of (3 dif- 
fers in different simulations from (3 = 1.0 10) to (3 = 1.5 
Q and may be non-universal for the objects of different 
mass scales . An attempt of analytical explanation of 
the results of numerical simulations has been performed 
m [HEI. The phase -space density profiles of DM halos 
are investigated in [l5| . 
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A central cusp in the Galactic halo and the smaller 
scale clumps result in the enhancement of DM annihila- 
tion rate and thus in the stronger signals in the form of 
gamma-rays, radio emission, positrons and antiprotons. 
The gamma-ray and radio signal from central cusp in 
the Galactic halo has been first discussed in Ila ll7i . Re- 
cently this problem was examined in [Tel ITfjl I20ll2ll |22| . 
The enhancement of DM annihilation rate due to the 
dumpiness of DM halo was first pointed out in pj . Neu- 
tralino annihilation in clumps can result in the very large 
diffuse gamma-ray flux |23| in the model of the clumpy 
DM by Gurevich et al. 07 Calculations of positron and 
antiproton production in the clumpy DM halo have been 
performed e. g. in 01 ( see a l so EH ES Eil ) ■ Recently 
the annihilation of DM in the clumps has been studied 

m mmmmoniiii. The s y ncnr ° tr ° n flux from ° m 

annihilation products in clumps in the presence of the 
Galactic magnetic field is considered in j 3 .' a] . Constraints 
on the DM dumpiness in the halos from heating of the 
disk galaxies is examined in |34Ll35|. 

The main purpose of this work is evaluation of the 
enhancement of the annihilation signal due to presence 
of the small clumps of DM in the Galactic halo. 

The small-scale self-gravitating Dark Matter Clumps, 
which will be referred to as DMCs or simply as clumps, 
may be formed in the early universe due to several mech- 
anisms. These DMCs may be formed (i) by the growth of 
the adiabatic or isothermal fluctuations (originated at in- 
flation) during the matter-dominated epoch; or (ii) from 
the density fluctuations in the models with topological 
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defects (cosmic strings and domain walls) or (iii) 
during the radiation dominating era from the nonlinear 
isothermal fluctuations (originated by phase transitions 
in early Universe) 13(1 or from the large amplitude adia- 
batic fluctuations [37|. 

In this paper we shall consider only the most conser- 
vative case of adiabatic fluctuations which enter the non- 
linear stage of evolution at the matter dominated epoch 
with the inflation-induced initial power-law power spec- 
trum. 

Small-scale clumps are formed only if the fluctuation 
amplitudes in the spectrum are large enough at the 
corresponding small scales. The inflation models pre- 
dict the power-law primeval fluctuation spectrum. If 
the power-law index n p > 1, DMCs are formed in a 
wide range of scales. During the universe expansion the 
small clumps are captured by the larger ones, and the 
larger clumps consist of the smaller ones and of contin- 
uously distributed DM. The convenient analytic formal- 
ism, which describes statistically this hierarchical cluster- 
ing, is the Press-Schechter theory [3^ and its extensions, 
in particular 'excursion set' formalism developed by Bond 
et al. ^39] (for the clear introduction see 40J ) . However, 
this theory does not include the important process of the 
tidal destruction of small clumps inside the bigger ones. 
We take into account this process in Sec. E] and obtain 
the mass function for the small-scale DMCs in the Galac- 
tic halo. In the case of the power-law spectrum only a 
small fraction of the captured clumps survives, but even 
this small fraction is enough to dominate the total anni- 
hilation rate in the Galactic halo. 

In the hierarchical theory of large-scale structure for- 
mation in the Universe the first formed objects have some 
minimal mass M m - ln . The value of this mass is determined 
by the spectrum of initial fluctuations and by the prop- 
erties of DM particles 0, ^| . This value is crucial for 
calculation of the DM annihilation rate. The estimates 
of M m i n existing in the literature for neutralino DM are 
substantially different, from M m ; n ~ 10 _12 M Q in to 
Mmin ~ (1CT 7 - 10- 6 )Mo in 43] ■ In Sec. ED we present 
our calculations and discuss the previous results. 

The DM annihilation rate crucially depends on the 
density profile p(r) oc r~ 13 of DM particles in a clump 
and on the distance R c where the density growth is cut 
off. This region is called the core. The radius of the 
core has been estimated in the literature in the different 
approximations. The estimation, R C J R ~ 5g , where <5 eq 
is the density fluctuation amplitude at the end of radia- 
tion dominated epoch, has been obtained in 0- It was 
found from behavior of the damped mode of nonlinear 
fluctuations. A black hole or baryonic core in the center 
of the DMC can strongly affect the density saturation at 
r -> for a very massive DMC H H Hfl EH . Calcu- 
lations 0, |2j| of the inward flux of DM particles into 
the dense central region of DMC also result, due to an- 
nihilation of DM particles, in a very small radius of the 
central core R c . The above mentioned process is essen- 
tial for the formation of the DMC core only in the case 



of almost perfectly spherically symmetrical clump. 

We shall estimate the radius of the core imposed by 
tidal interaction, which gives the largest R c among those 
known in the literature. In the spherically symmetric 
self-gravitating clump on the stage of its formation the 
non-dissipative DM particles are moving nearly radially 
in the oscillation regime. The presence of a non-spherical 
(tidal) external gravitational field causes the deflection 
of particle trajectories in the clump from the radial ones. 
This process prevents the development of the central sin- 
gularity in a clump and results in the core formation. 

During the radiation dominated epoch the small fluc- 
tuations 5 = Sp/p <C 1 grow very slowly. At the matter 
dominated stage these fluctuations start to grow fast 01 
in the regime 6 cx i 2 / 3 . Fluctuations get detached from 
the general cosmological expansion and start contracting 
after reaching the nonlinear value, 5 > 1. These nonlin- 
ear fluctuations form finally the DM clumps. The ana- 
lytic studies of nonlinear evolution of fluctuations have 
been performed by many authors. One of the most de- 
tailed analytical approach was developed by Gurevich 
and Zybin [9|. In this formalism at the certain moment 
of gravitational contraction the density singularity forms 
in the center of a nonlinear fluctuation. From this sin- 
gularity point the density caustics (i. e. the boundary of 
regions with a different number of streams) expand out- 
ward. The secondary caustics appear inside the primary 
ones and their number increases fast with time. This 
multi-stream instability has been discovered and studied 
in detail in |9|. It was demonstrated that the stationary 
universal density profile p(r) oc r^ 13 with (3 = 1.7 — 1.9 
is formed as a result of streams mixing. The maximum 
density of DM particles in a clump is reached at the cen- 
ter. 

In our consideration of the clump formation we shall 
follow for convenience the theoretical scenario of Gure- 
vich and Zybin 0. However, the effects of tidal interac- 
tion, which is the main result of our work, are valid for 
a much broader class of scenarios. 

The processes described above are valid for all DM 
particles which can be considered as non-dissipative. 
The signal production depends on the annihilation cross- 
section and thus on the nature of DM particles. However 
our strategy is to calculate the enhancement of the sig- 
nal due to the halo dumpiness in comparison with an 
isotropic unclumped distribution of DM. As a guide we 
shall take neutralino as DM particle, but essentially our 
results for enhancement of the annihilation signal are rel- 
evant for a wide class of other DM particles. 

We perform our calculations for cosmological model 
with the matter density 17 m = 0.3 and the cosmologi- 
cal term Oa = 1 — Q m ~ 0.7. The presence of A term 
influences only the value of p cq and does not affect the 
formation of low mass DMCs. This is because the A term 
contributes negligibly to the total cosmological density at 
time scales when the low-mass DMCs formation occurs. 
We shall use index 'eq' for the values at the moment of 
equality (i. e. transition from the radiation dominated 
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to matter dominated epoch). We shall use the Hubble 
constant 70 km s _1 Mpc -1 . 



II. ENHANCEMENT OF ANNIHILATION 
SIGNAL DUE TO CLUMPS 

Let us consider a DM clump with the internal density 
profile Pmt(f) and a total mass M = J 4irr 2 p lnt {r)dr. An 
annihilation rate in a single clump is given by 



iV c i = 47r r 2 drp 2 lt (r)m 2 (a l!inn v) = 



4tt 



R 3 



S, 



(1) 

where m x is mass of a DM particle (being not necessarily 
neutralino), v is relative velocity of two DM particles 
at the collision, er ann is annihilation cross-section and R 
is virial radius of a clump. Function S is determined 
by Eq. Q and depends on DM distribution in a clump, 
in particular S = 1 for the simplest case of an uniform 
clump, when pi n t(r) = const at r < R and /Oj n t(f") = at 
r>R. 

An expansion of (<7 an n^) over the relative velocity v of 
two DM particles has a form 



(cr a 



cv 



(2) 



where a has a contribution of s-wave amplitude only, and 
b — from both s- and p- waves. Since v is very small, 
(<7ann^) can be put out of the integral in Eq. (Q. 

We shall use the following parametrization of the den- 
sity profile in a clump 



PintO) 




r<R r 



-0 



R~J ' R c <r<R > 
r> R. 



(3) 



Using p- m t(f) from Eq. (J3J) it is easy to calculate S from 
Eq. 0) as 



S(X C) 0): 



(3-/?) 2 {2/3 



3(2/3-3) 



• -,-: ! 1 



(4) 



where x c = R c /R. Another approach to the parametriza- 
tion of clump structure has been used in [23 ■ 

There is distribution of clumps in the Galactic halo 
at least over three parameters, mass M, radius R, and 
distance from the Galactic Center I : n c \(M, R,l). This 
distribution can also depend on the parameters which 
describe the internal structure of the clumps, j3 and x c , 
from Eq. 1(3}. We shall discuss this dependence in Sec. II VI 
In particular it will be demonstrated that x c — x c (M, R), 
while j3 is the universal constant. Thus the differential 
number density of DMCs in the halo can be written as: 



dN r 



(l,M, R)d 3 ldMdR, 



(5) 



The observed signal at the position of the Earth from 
DM particle annihilation in the clumps is proportional 
to the quantity 

7T r maI (C) M max _R max 

7 cl = i- | d CsinC / ^ JdM jdR 

M min R mi „ 



x n cl (l(C,r),M,R)N cl (M,R), 



(6) 



where r is distance from the Sun (Earth) to a clump and 
C is angle between the line of observation and the direc- 
tion to the Galactic center. The distance I between a 
clump and the Galactic center can be given in terms of 
r, tq (distance from the Sun to the Galactic center) and 
£ as l((,r) — (r 2 + Tq — 2rr Q cosC) 1 ^ 2 , and the maxi- 
mum distance from the Sun to the outer halo border in 
the direction of (, r max (C) = (R 2 H — sin 2 C) 1//2 , where 
Rh ~ 100 kpc is the Galactic halo virial radius and 
r Q = 8.5 kpc is the distance from the Sun to the Galactic 
center. 

Additional annihilation signal is given by unclumpy 
DM in the halo with homogeneous {i.e. smoothly spread) 
density pdm(0j where / is a distance to the Galactic Cen- 
ter. 



dC, sin £ 



c(0 

drp 2 



DM 



(2(C,r))/m 2 . (7) 



The enhancement r\ of the signal due to a presence of 
clumps is given by 



V = 



Id + Iho 



Iho 



(8) 



This quantity describes the global enhancement of the 
annihilation signal observed at the Earth (e. g. the flux 
of radio, gamma, and neutrino radiations) as compared 
with usual calculations from annihilation of unclumpy 
DM. From Eqs. ©, and © one can see that en- 
hancement T) does not depend on the properties of DM 
particles, in particular on the annihilation cross-section, 
and is fully determined by the parameters of DM dumpi- 
ness. The further exact calculations in this paper will be 
performed using Eqs. JSJ-JSJ, but now we shall turn to 
the approximate expression for n. 

We shall accept now the simplifying assumptions. We 
assume that space density of clumps in the halo, n c \(l) 
is proportional to the unclumpy DM density, pum(1)' 
n c\(l) = £pdmQ)/M with (< 1. This assumption holds 
with a good accuracy for the small-scale clumps. In con- 
trast, the distribution of large-scale clumps obtained in 
the numerical simulations |45j is rather different from the 
density distribution of the small clumps, especially in a 
central part of the halo because of the tidal disruption of 
clumps there. However, the clump signal is determined 
mostly by clumps of the minimal mass. We neglect here 
the distribution of clumps over M and R. Instead we 
shall use a mean density of DM particles inside a clump 
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Pint = 3M/47ri? 3 . Finally, we shall introduce for conve- 
nience the effective density of DM particles in the halo 
defined as 



Pdm 



(9) 



»w(0 

dCsinC J drpn M (l((,r)) 
o o 

r max (C) 

dCsinC J drp DM {l((,r)) 
') o 



As a result, we obtain for an enhancement factor the 
convenient estimate 



Pdm 



(10) 



where £ is a fraction of DM in the form of clumps (see 
above) and S(x c ,f3) is given by Eq. J3J. For typical pa- 
rameters (see details in the following sections) n p ~ 1, 
/3 ~ 1.8, x c ~ 0.05, S(x c ,/3) ~ 5, pjj ~ Pdm(^ ) 



0.3 GeV cm- 



Pint 



2 x 10~ 22 g cm" 3 , £ - 0.001, the 



numerical estimate 7/ ~ 3 follows from Eq. 110J1 . 



III. CLUMPS OF MINIMAL MASSES 

The number of clumps in the halo increases at small 
clump masses M, and the signal from clumps I c \ crucially 
depends on M m ; n in the clump distribution as Eq. © 
shows. The value of M m i n is determined by a leakage of 
DM particles from the overdense fluctuations in the early 
universe. We shall describe first this process qualitatively 
and then present numerical calculations. 

CDM particles at high temperature T > Tf ~ 0.05m x 
are in the thermodynamical (chemical) equilibrium with 
cosmic plasma, when their number density is determined 
by temperature. After freezing at t > tf and T < Tf, the 
DM particles remain for some time in kinetic equilibrium 
with plasma, when the temperature of CDM particles T x 
is equal to temperature of plasma T, but number density 
n x is determined by freezing concentration and expan- 
sion of the universe. At this stage the CDM particles are 
not perfectly coupled to the cosmic plasma. Collisions 
between a CDM particle and fast particles of ambient 
plasma result in exchange of momenta and a CDM parti- 
cle diffuses in the space. Due to diffusion the DM parti- 
cles leak from the small-scale fluctuations and thus their 
distribution obtain a cutoff at the minimal mass Mjj. 

When the energy relaxation time for DM particles r rc i 
becomes larger than the Hubble time -ff -1 (t), the DM 
particles get out of the kinetic equilibrium. This condi- 
tions determines the time of kinetic decoupling t&. At 
t > td the CDM matter particles are moving in the free 
streaming regime and all fluctuations on the scale of 



a(t ) 



t 



v(f)df 
«(*') 



(11) 



and smaller are washed away [here a(t) is the scaling 
factor of expanding universe and v(t) is velocity of a DM 
particle at epoch t] . The corresponding minimal mass at 
epoch to, 



4tt 

T 



p x (t )X%, 



(12) 



is much larger than M]j. Numerical calculations below 
(for neutralino) show that M/j is close to M m ; n from 
and Mf s to M m j n from |43| . 

The calculation of the minimal mass Mmin in the mass 
spectrum of fluctuations is obviously model dependent. 
As the DM particle we shall consider the neutralino \i 
for which we take the pure bino state (x = B). As calcu- 
lations below show the temperature of kinetic decoupling 
for a reasonable range of parameters is Td ~ 100 MeV, 
and thus we can consider cosmic plasma consisting of 
relativistic electrons, positrons, neutrinos and photons 
in thermal equilibrium. 

The cross-sections for scattering of bino off left (right) 
electron and left neutrino are given in the Appendix A. 
A cross-section for v\ scattering is given by Eq. ilA5(l and 
for e\ scattering it is 17 times larger, if to assume equal 
masses of selectrons and sneutrinos (we shall use m for 
the both left and right selectron and sneutrino masses, 
and M 2 = m 2 - to 2 .). 

First of all we shall calculate the moment td and tem- 
perature Td of kinetic decoupling of neutralino, using con- 
dition 



— ^H(t), 



(13) 



where H (t) = 1/ (2t) is the Hubble constant and r re i(T) is 
the energy relaxation time for neutralino at temperature 
of electron-neutrino gas T. The relaxation time t ic \ is de- 
termined by collisions of neutralino with fermions vl, &l 
and eR. Neutralino can be considered as particle at rest 
because the rest system coincides with the center-of-mass 
system with the accuracy of order y/T/m x . Let 5p is the 
neutralino momentum obtained in one scattering: 



(Sp) 2 = 2w 2 (l - cos 6»), 



(14) 



where u> and 9 is neutrino energy and scattering angle, 
respectively. 

Let us introduce the number density of relativistic 
fermions with one polarization and with energy lo: 



n (w) 



1 uj 1 



2tt 2 e"/ T + 1 ' 
Then for the energy relaxation time r ro i we have 



(15) 



1 1 dEk 40 [dnfdvHoUlg) (s P ) 2 



T rc l Ek dt 2£fcTO a 



fLX 



(16) 

where Ek — (3/2)T is a mean kinetic energy of neu- 
tralino, and (da e \/ dVt) f LX is given by Eq. i|A5|l . The 
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number 40 in Eq. (|16fl is obtained by counting of de- 
grees of freedom: three neutrinos and antincutrinos (or 
vl in case of Majorana neutrinos) give 6, and e c L give 
2 and two right (singlet) states for electron and positron 
gives 34, because their cross-sections are 17 times larger. 
After integration in Eq. (|16|) we obtain 



1 = 40r(7)a e 2 m . T & 
r rc i 9ttcos 4 6> w M 4 m. 



(17) 



Using Eq. I|13|) and connection between age and temper- 
ature of the universe 



t = 



2.42 



T 



1 MeV 



(18) 



where g* is number of degrees of freedom, we obtain nu- 
merically: 



t d = 3.5x10" 



and 



m x \-V2/ M \ /g»\- 3 / 4 



V100 GcVJ 



1 TcV 



(19) 




MeV. 



(20) 

We shall present in this section the calculations made 
in physically transparent way, considering the diffusion 
leaking of neutralinos at the stage of kinetic equilibrium 
and free streaming when neutralinos get out of kinetic 
equilibrium. In Appendix B we shall study both stages 
together in the formalism of kinetic equation, as it has 
been done in [42|. Though our methods are not iden- 
tical, their comparison implies that the absence of free 
streaming is responsible for the contradiction with differ- 
ent Mmin discussed above. The independent approach in 
Appendix B confirms the results obtained below. 



A. Diffusion cutoff of the mass spectrum 

We can come over now to the calculation of Md, the 
minimal mass in the fluctuation spectrum caused by dif- 
fusion of neutralinos out of an overdense fluctuation. We 
calculate the diffusion coefficient using the method given 
in |46l | (§12). Consider a neutralino moving with a non- 
relativistic velocity v. In the rest system of this parti- 
cle the momentum distribution of relativistic fermions is 
anisotropic: 



n(p) d 3 p - 



dfl a p 2 dp 



1 



(2tt) 3 e p(l+« cos a)/T + 1 - 



(21) 



where a is the angle between the directions of v and 
momentum of incoming fermion. 

The momentum transfer in a single scattering equals 
to p{l — cos 9) after averaging over the azimuthal angles. 



A corresponding friction force experienced by the neu- 
tralino is 



f r = 40 Jdilg I d 3 pn(p) 
= -B^v, 



d(T c _ 

dK 



p(l — cos ( 



./'/ X 



(22) 



where B is a particle mobility and factor 40 takes into 
account scattering on all fermions as in Eq. i|17|) . Then 
the diffusion coefficient is 



D = BT = 



3?rcos 4 6» W M 4 



40r(6)a| m T5 
Diffusion equation in the comoving system has a form 
8S(x, t) D{t) 



(23) 



dt a 2 (t) 



(24) 



Diffusion coefficient D(t) is time-dependent because of 
T(t). Solution of Eq. 12411 for the Fourier component is 

6 n (t) = 8 % {t } ) exp [-k'Cgl^M* (i 5 ' 2 - t 5 / 2 ) } , (25) 

where C = const. The factor Cgl /4 M 4 (t 5 / 2 - t 5 / 2 ) in 
front of k 2 in Eq. (|25|l is the diffusion length squared 
Aq(<) / 'a 2 ; (t) in the comoving coordinates. This value de- 
termines the minimal mass in the density perturbation 
spectrum due to diffusion of neutralinos from a fluctua- 
tion: 



Mr. 



4?r 



M^)A^)=4.3xl0- 13 ( 



m x \-l5/8 



100 GeV 



/ ~ \ -3/2 
I M \ /Q*\- 15 /16 

-b (is) < 26 > 



The functional dependence of Eq. 125fl and numerical 
value of l|2()l) obtained in diffusion approximation coin- 
cide with the corresponding results obtained by different 
method in 421. 



B. Free streaming cutoff of the mass spectrum 

We shall consider now the free streaming cutoff of the 
mass spectrum qualitatively described in the beginning 
of this section. We have given there an estimate of the 
minimal mass due to free streaming. In the accurate 
calculations below we shall take into account the angular 
and velocity distribution of leaking neutralinos, and exact 
dependence of a(t) at age ~ £ eq , which affect the value of 
M fs . 

After the moment of kinetic decoupling td, neutrali- 
nos move freely in the expanding universe background, 
a(t)dx — v{t)dt, where x is comoving particle coordi- 
nates. Coordinates x at some moment t are determined 
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by initial coordinates q and velocity v*d at the moment of 
kinetic decoupling td'. 



resulting M m - ln at to can be easily calculated using a(t) 
as solution of the Friedman equation: 



x = f(q, Vd, t) = q+ I " ( *i/f = 9 + 9(t)v d , (27) 



a(t>) 



where 



g(t) = a(t d ) 



df 



a 2 {t')' 



(28) 



v(t) = Vd<x(td) I 'a(t) for nonrelativistic particle. Now for 
the neutralino number density at the point x we have 



n(x,t) = Jd 3 v d (f>{vd)^2n(q*,t d 



£1 

Dq 



(29) 



d 3 v d (j)(vd) / d 3 q n(q,t d )S i3) (x- f(q,v d ,t)), 



where <5^ 3 ^ is the Dirac delta- function, Df /Dq is the Ja- 
cobian and 4>{v d ) is neutralino velocity distribution func- 
tion at the moment td- Summation in H29|) goes over 
all roots <f» of the equation x — f(q,v d ,t) from l(77)l. 
This sum in fact has only one term because the function 
f(q, Vd, t) in f?7| is a single- valued one. 

From l|27|) and (|29[1 we find the Fourier component 



(30) 



Assuming velocity distribution at the moment of decou- 
pling 4>{v*d) to be Maxwellian 



= ( m x A 



3/2 



[Vd) = {^Fd) exp t 2T d r 



we obtain 



n k (*) = n fc(^) e " 



(31) 



(32) 



i.e. up to the moment i all perturbations are washed out 
by free streaming inside the physical length-scale 



A fe (i) = a(t)g(t) ( ^ 



1/2 



(33) 



This length-scale corresponds to the clump of the mini- 
mal mass 



47T 



(34) 



where = /0 C q»cq/a 3 (i)- The time dependence of 

Mf s (t) is regulated by a(i). At the radiation dominated 
epoch, Mf s (t) grows logarithmically with time. This 
growth is saturated at the matter dominated epoch. The 



t dcqTJ* 



I (l 

3 V, 7 ?* 



(35) 



In these equations i]^ 2 = 2irGp cq a 2 q /3, a cq is the value 
of scaling factor at the moment i eq , 

Pcq = Po(l + z cq ) 3 = 1.1 x 1(T 19 ( AV cm -3 ; 

(36) 

1 + z oq = 2.35 x 10 4 ft m /i 2 and p = 1.9 x 
10 -29 r2 m /i 2 g cm -3 . The presence of cosmological con- 
stant A affects only the value p eq and does not influence 
the evolution M{ s (t) because the contribution of A to the 
total cosmological density is negligible at small t. Putting 
(|35|l into Q28p. we find after integration 



_l/4 1/4 y.3/2 / rn \ 3/2 



2 19/4 3 l/4 £3/4 v m 



24 



TTGp eq t d 



Using Eqs. I|19|) and (|20[1 we obtain numerically 



(37) 



M„ 



1.5 x 10" 



( m x \ 
V100 GcVJ 



-15/8 



71/ 



-3/2 



1 TcV 



ff»\- 15 / 16 /A 



(I) 



83' M0 ' 



(38) 



where A* is the logarithm from Eq. (I37|l . 

Our calculations agree well with |43| as far as the most 
important quantity Td is concerned (the scattering cross- 
section is involved only there). The calculation of M m ; n 
from Td in our case involves non-radial propagation of 
neutralinos in a fluctuation and their distribution over 
velocities, Eqs. (|27H - <|32[) . We include also the accurate 
time dependence of the scaling factor ait). The calcu- 
lations in 43] follow the semi-quantitative scheme de- 
scribed in the beginning of Sec. II I II At this stage of 
calculations we have a difference described by factor 7. 

In conclusion, in this section we have considered two 
processes of washing out the cosmological density pertur- 
bations in neutralino gas. The first process is the neu- 
tralino diffusion due to scattering off neutrinos, electrons 
and positrons. This process is effective until neutrali- 
nos are in the kinetic equilibrium with the cosmological 
plasma. Up to the moment of decoupling td all pertur- 
bations with mass M < M D ^ 10~ 13 - 1O~ 12 M are 
washed out. The second process is neutralino free steam- 
ing. Starting later, at t > t d , it washes out the larger 
perturbations with M < Mf s and determines M m j n in 
the clump mass distribution at present epoch, as given 
by Eq. (EEJ. 
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IV. CORE OF A DARK MATTER CLUMP 

In this section we shall consider smearing of the sin- 
gular density profile in a clump due to tidal forces and 
calculate the radius R c of the produced core. 

Clumps, as well as galaxies, are originated near the 
maxima in cosmological density perturbations S(r) = 
(p(r) — p) I p. At the matter dominated stage the density 
perturbations grow as 5 oc t 2 ^ 3 . In the nonlinear regime, 
5 > 1, the multi-flux instability develops in a clump and 
singular density profile is formed 9]. If velocity field in 
the central part of the clump is disturbed and becomes 
weakly nonradial, the flow is overturned, singularity does 
not form and density profile is smoothed. In the core 
radius is estimated as x c ~ S 3 q <C 1 from consideration 
of the perturbation of the velocity field due to damped 
mode of the cosmological density perturbations. Here S oq 
is an initial density fluctuation value at the end of radi- 
ation dominated epoch. In the core is produced for 
spherically symmetric clump by inverse flow caused by 
annihilation of DM particles. We shall show here that 
these phenomena are not the main effects and that much 
stronger disturbance of the velocity field in the core is 
produced by tidal forces. These forces originate due to 
non-sphericity of the considered perturbation and pres- 
ence of other fluctuations nearby, including a fluctuation 
of larger scale in which the considered fluctuation can be 
submerged. 

The core formation in a fluctuation begins at the lin- 
ear stage of evolution and continues at the beginning of 
non-linear stage. The tidal forces diminishes with time 
t (see Eq. (I53fl below), while duration of this phase is 
proportional to t. Once the core is produced it is not de- 
stroyed in the evolution followed. The stage of the core 
formation continues approximately from t eq to the time 
of maximal expansion t s and a little above, when a clump 
is detached from expansion of universe and evolves in the 
non-linear regime. Soon after this period, a clump enters 
the hierarchical stage of evolution, when the tidal forces 
can destroy it, but surviving clumps retain their cores. 

Let us expand the gravitational potential in the series 
near the maximum of the density fluctuation taken as 
r = at arbitrary time t during the linear growth of 
density perturbations: 



0(f,i) 
where 



dj)_ 

Q r l 



1 



1 



r l + - S H | %rV+-T ( y| rV + 



d 2 (j){r) 
dr l dr^ 



Si 



:<$>llSi 



(39) 



(40) 



The first term of series in Eq. I|39[) does not influence 
the particle motion. The second term is zero as a condi- 
tion of maximum density. The third term describes the 
spherically symmetric part of the potential (including the 
potential of the homogeneous background) and also the 
perturbation potential. It governs the radial motion of 



the particles. According to the Poisson equation one has 

S w = A0(r) = 47rGp(I + 5(f)). (41) 

Finally, the fourth term, which contains the traceless ma- 
trix Ty, describes the tidal forces. They disturb the ra- 
dial motion of the particles and result in production of 
the core. 

We shall start with definitions and notation. We as- 
sume that density perturbations S(r) are Gaussian with 
a power spectrum P(k): 

(Stfz) = (2n) 3 P(k)6^ (k - k'), ^ = J S(r)^d 3 r, 

(42) 

where 5^ (k — k') is the Dirac delta-function and angle 
brackets corresponds to ensemble averaging. The power 
spectrum P cq (k) at the time t cq is connected with the 
primordial power spectrum P p (k) (at the epochs before 
the horizon crossing) by relation P cq (k) = P p (k)T 2 (k), 
where T(k) is the transfer function for cold dark matter 
(see e. g. |47|'). 

From Eq. (|4I|1 it follows that the power spectrum 
P$(fc) of potential perturbations is related to P{k) as 

P^(k) = {4ir) 2 G 2 p 2 k~ 4 P(k). (43) 

Let us introduce the moments of spectrum P{k) 

oo 

I 



r (i) 



2vr 2 



k 2 dkP(k)k 2j , 



(44) 



and the similar moments s 2 ^ for the perturbation field 
of the gravitational potential. Calculating the divergent 
moments for given mass M we assume smoothing proce- 
dure of [13. 

From Eq. (g3J it follows that 



(4n) 2 G 2 p 2 a 2 



0-2) 



(45) 



for j > 2. Let us define Cij = d 2 S \r) / dr % <9r? . Then 
according to |47]|. its mean value over the ensemble is 



( 2 ) 

(djCki) = -^-(SijSki + SikSji + SuSjk), 



which results in 



(TijTji) = 3 S (2) = g( 47r ) 2G Ve r ( ! o)> 



(46) 

(47) 
Let us 



(in the following we shall use a notation a = ct( ) ) 
introduce the important physical quantity ^, the height of 
the peak density in units of dispersion (the peak- height): 



v = 5 cq /a cq (M), 



(48) 



where a cq (M) = a(t cq , M). 

After introduction of these quantities we shall come 
over to calculation of DM particle velocities and core ra- 
dius. The velocity v(t) is given by sum of radial velocity 
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u ra d and an additional velocity Vtid, which will be con- 
sidered as a small correction caused by tidal forces. The 
radial velocity will be calculated without tidal interaction 
taking into account from the equation 



dv T 



.•I 



dt 



= -grader), 



(49) 



where spherically symmetric potential <j)(r) is given by 
the third term in r.h.s of Eq. Ij39|l . This equation deter- 
mines the radial motion of the particle, and its solution 
is given in the parametric form as |44j |: 



1 



sin 26 



5<5 C 



3/2 



t-t a 
tea 



(50) 



The moment of maximum clump expansion t s and the 
distance r — r, at this moment are 



if 

tec 



3tt 
T 



56, 



cq 



-3/2 



Tec 



5<5 C 



(51) 



where <5 eq is the initial fluctuation value (at t = t cq ). 

Tidal forces give rise to the additional velocity Vtid- Its 
evolution is described by equation 



dvtic 



± = - Tl] (ty 



dt 

where in the linear approximation 

t 



Tij(t) — Tij(t eq ) 



cq 



-4/3 



(52) 



(53) 



The linear approximation for tidal forces is justified be- 
cause they are generated by the large-scale perturbations 
which become nonlinear later than the small-scale pertur- 
bation under consideration. 

Now we find the value Vtid at the moment when the 
density nonlinearity sets in, S ~ 1, or more exactly at the 
moment of a maximum expansion t s . After integration 
of with the help of Eq. J3UJ) we obtain 



vtid,i(* s ) = 18 1/3 f -^J t cq f(S cq )T t] (t cq y(t s ) 



where the function 

tt/2 
(55 cq /3)!/2 



(54) 



■ sin 2(f) 



-4/3 



sin 4 (55) 



is plotted in the Fig. ^ We may use approximately 
/(<5 oq ) — 1 for the values of <5 e q at interest (asymptot- 
ically f(6 eq ) -> 1.32 at 5 cq -> 0). 

To find the core radius R c of the clump we shall use 
the method similar to that in 9] . Since rot v t id = and 
divu ti( i = 0, the tensor has the following diagonal 
form in the coordinate system connected with the main 
axes. 



log S eq 

FIG. 1: Function /(£«,) defined by Eq. The ratio of the 

core radius to the clump radius in typical case ~ 0.3v~ 2 f(6 cq ). 



T' = 



-2t 



(56) 



The value r from Eq. I|56l) is connected with the core ra- 
dius R c due to energy relation AE ~ AV, where AE ~ 
J d 3 rpi n t(t s )vt id (t s ) /2 is the work performed by tidal 
forces, and AV ~ GMMJR, where M c ~ 4irRlp c /3, 
is potential energy of the core. It gives for the relative 
core radius 



2 2/3 3 10/3 T 2 t 2 Si 



cq^cq 2 



(57) 



From invariance of the matrix trace relative to the change 
of coordinates one has 

r 2 = I*.2£./6 = TijTji/6. (58) 

Since the correlator 

(TijS) = (($„■ - i$„%)($ H (4 7 rGp)- 1 - 1)) - 0, (59) 

the quantities S and TL are statistically independent and 
we may average them independently. Averaging (r 2 ) over 
the tidal force field Ty with help of Eqs. (03) and JSBJ), 
we obtain the main result of this section for the relative 
core radius: 



Rc 
R 



^5/3313/3 



.., Gp cq t 2 cq v 2 f 2 (S cq ) 



0.3^ 2 / 2 (<5 e q), 



(60) 



where v is given by Eq. (j48(l . 

The fluctuations with v ~ 0.5 — 0.6 have x c ~ 1, i.e. 
these fluctuations are practically destroyed by tidal in- 
teractions. Most of galactic clumps are formed from 
v ~ 1 peaks. As it will be demonstrated in Sec. IVII 
those clumps that survive till now are characterized by 
v ~ 1.6, but the main contribution to the annihilation 
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signal is given by the clumps with v ~ 2.5 for which 
x c ~ 0.05. 

In an alternative approach for calculation of the core 
radius one may define R c as the minimum deviation of 
a typical particle trajectory from a center of the clump. 
After clump virialization a particle at the average dis- 
tance R/2 from the center has an angular momentum 
mi>tid-R/2. At t > t s the tidal force is already small and 
angular momentum is approximately concerved. Then 
at the minimum distance from the core, R c , one has 
R C V ~ v tid R/2, where V = {2GM/R) 1 ' 2 is the velocity 
of typical particle in the center of the virialized clump. 
Calculating w t id from Eq. ijMj) and using Eq. (|T7jl one 
obtains 



X c ~ O.IS^/OW, 



(61) 



which numerically is very close to Eq. 1|60|) for typical 
values of v ~ 1 — 3. 

The core radius, given by the Eq. (|60|) . is much larger 
than x c ~ 5g q obtained in |j| and x c from .23j. The core 
radius found in ;23] is valid only in a perfectly symmetric 
case. 

Several remarks are in order. 

Tidal forces prevent the appearance of singularity dur- 
ing evolution of the clump, but if such singularity has 
somehow appeared, tidal interactions cannot destroy it. 

In the calculations above we operated with average 
tidal force, described by Eq. Q47[l. In reality this force 
fluctuates due to positions of surrounding fluctuations 
which can overlap with the considered one or to be far 
away from it. As a result some clumps can be destroyed 
and those survived have different core radii R c . This 
effect increases the annihilation signal. 

We assume that DM distribution within the core is 
much flatter than r~@. Between R c and the beginning 
of the asymptotic regime /?i n t oc r~@ there is a transi- 
tion zone. During the hierarchical evolution (see next 
section) this zone expands due to tidal interaction in the 
hierarchical structures. However, this interaction cannot 
destroy the existing core. 

The above calculations are also valid for the formation 
of galaxies. It is useful to compare Eq. (|60|) with obser- 
vations of galaxies. In ;47j the number of peaks in the 
Gaussian random field is confronted with the observed 
density of spiral and elliptical galaxies. It was found 
that these galaxies are formed mostly from peaks with 
v ~ 3. According to Eq. lIBTijl for these peaks x c ~ 0.033. 
The rotation curves in the central part of dwarf and Low 
Surface Brightness galaxies are measured |4Sl 149^ and 
constant-density cores were revealed. Existing observa- 
tions do not contradict the presence of extended uniform 
core with radius R c = x c R ~ 3.3 kpc in spiral and ellip- 
tical galaxies. However, at these distances the baryonic 
matter dominates, which makes the observation of DM 
core even more difficult. The extreme value R c = R 
corresponds to t-Wn « 0.55/. The Gaussian peaks with 
these v are completely washed out by tidal forces and 
do not produce the gravitationally bound objects. The 



intermediate case v ~ 1 corresponds to most numerous 
dwarf and irregular galaxies with the pronounced over- 
density in the central core. 



V. TIDAL DESTRUCTION OF CLUMPS IN 
THE HIERARCHICAL MODEL 

In this section we shall study the destruction of clumps 
by the tidal interaction which occurs at the formation 
of the hierarchical structures, but long time before the 
galaxy formation. This interaction arises when two 
clumps pass near each other and when a clump moves 
in the external gravitational field of the bigger object 
(host) to which this clump belongs. In both cases a 
clump is exited by the external gravitational field, i. e. 
its constituent particles obtain additional velocities in the 
c. m. system. The clump is destroyed if its internal en- 
ergy increase AE exceeds the corresponding total energy 
\E\ ~ GM 2 /2R. In Sec. lVBl we shall calculate the rate of 
excitation energy production due to both aforementioned 
processes. Respectively in Sec. IV Cl we shall calculate the 
survival probability for a clump in the hierarchical model, 
in which the smaller clumps are embedded in the bigger 
one, and the latter enters into a more bigger clump etc. 
But first we shall describe the necessary generalities and 
definitions. 



A. Generalities and definitions 

The formation of DM objects with a fixed mass M 
at the linear stage is distributed over formation epochs 
t f . In the spherical model of the Press-Schechter theory 
|38l l40| the density perturbation at the epoch of object 
formation is equal to S c = 3(12tt) 2 / 3 /20 ~ 1.686: 

S(M,t f )=S c (62) 

Eq. (|62|l gives the formation criterion for DM objects. 
The formation criterion alone does not determine tf for 
a given mass M, because 6(M) has the Gaussian distribu- 
tion. The formation time t f can be fixed by an additional 
condition, e. g. v = 1 (sec Eq. (0HJO- The DM objects 
which satisfy the formation criterion (|62l) and v = 1, 
or equivalently a(M,tf) = S C) will be referred to as the 
typical objects. For a given mass M they are character- 
ized by a fixed epoch of formation tf. In some parts of 
our consideration we shall simplify the problem, assum- 
ing that clumps are typical instead of taking into account 
their distribution over tf. 

We confine ourselves here only to the power-law spec- 
trum of fluctuations P oq (k) oc k n , in which case 



CTcq (M) cx M-(™+3)/6. 



(63) 



The effective power-law index n obtained from the ex- 
pression above can be given as 



-3-6 



dln<r cq (M) 
dlnM ' 



(64) 



10 



In the case of an arbitrary P(k) spectrum, the a eq (M) 
has also an arbitrary dependence on M. Eq. l|64l) defines 
n e ff for a some mass interval AM. In realistic cases (see 
e. g. Eq. iJOU), Eq. (|S4"|) is a good approximation because 
a power-law index n depends only weakly on M. 

We shall introduce for convenience the characteristic 
values of typical objects labeled by index 'A' and de- 
scribed by the condition a(t\,M\) — S c at some fixed 
moment t\ or redshift z\. We can choose these quanti- 
ties arbitrarily (because they will not enter into the final 
results), but satisfying the condition za S> (p,A./ft m ) — 1, 
e. g. taking za ~ 5 — 10. The convenience of these 
normalizing values is that at t < t\ the A-term can be 
neglected, while small-scale object formation occurs at 
these epochs. Let us introduce also the dimensionless 
mass m as 

m = M/M\. 
Using the formation criterion (|62J) we obtain 



S(M, t 



f) 



Sea 



1 



1 



= S r . 



Z f 



(65) 



because the growth- factor for the rising mode D g (t) oc 
(1 + z) cx t 2 l z in the standard cosmological model at the 
matter dominated epoch. For a single clump with mass 
M obeying S cq = va eq (M) with an arbitrary v one has 



Pint = Kp(Zf) = np cq 



l + z 



eq 



np c 



1 51 



and 



R 



( 3M 



1/3 



i?ATO 



(n+5)/6) 



(66) 



(67) 



where M and R are respectively mass and radius of the 
clump, k = 18tt 2 ~ 178 and R A = (3M/i7rKp(t A )y/ s . 
The value of k describes the clump density increasing 
during contraction and it can be found from Eq. 15 01) 
(see also [iojl. 

A number density of unconfined (free) clumps (i.e. of 
those not belonging to the more massive objects) is given 
by the Press-Schechter formulae jijj 



1/2, 



M D g (t)a 2 q dM 



x exp 



b PS (t,M)dM 
-S 2 



2D g (t) 2 a 2 q 



dM, 



(68) 



where the growth- factor D g (t) is normalized as D g {t eq ) = 
1. Let us consider the Press-Schechter distribution for the 
clumps hosted by the larger clump at the moment tt of 
its formation. Taking into account the density increase 
factor k we obtain 



ipps(t f i m ) dm = K<pps(tf,m)dm. 



(69) 



As it will be demonstrated in this section the clumps 
are destroyed by tidal interaction and each clump has a 
small surviving probability £ < 0.01. A survived clump 
during its lifetime is surrounded by other clumps with 
distribution given by Eq. (|69() . When a host clump is de- 
stroyed, the survived clump finds itself hosted by larger 
clump with the small clump distribution inside given by 
the same Eq. (|69f) but with larger t / . Since the tidal de- 
struction is most effective at small distances, one should 
always consider the smallest possible host clump from the 
hierarchical structure, and the distribution of the small 
clumps around one at the consideration is always given 
by Eq. (|69l) . The characteristic time is the time of for- 
mation tf of the smallest host clump, though time of 
destruction can be much larger than t f . 

The total energy (kinetic and potential) of a clump is 
given by 



\E\ 



3-/3 GM 2 
2(5 - 2/3) R 



(70) 



B. The rate of internal energy growth 

Consider a host clump with mass Mh, radius Rh, and 
with the small clumps inside having distribution (|69|l and 
moving in the common gravitational potential with a ve- 
locity dispersion ~ Vu — GMh/Rh- Interacting tidally 
with its neighbors, a small clump increases its internal 
energy. We calculate first the rate of internal energy in- 
crease due to these interactions. The mass of the consid- 
ered clump is M = toMa, it is characterized by an arbi- 
trary v and its interaction with a target clump occurs at 
the impact parameter I. A target clump is characterized 
by mass M' , radius R', radius of the core R' c = x c R with 
x c ~ 0.1 and by universal density distribution 0. 

An increase of internal energy of a clump M during 
one fly-by in the momentum approximation [50| is given 
by: 



1 



r, ^2 



AE=- d 6 rp{r){v x -v x ) 



(71) 



where v x is the velocity increase of DM particle in di- 
rection of axis x, and v x is the same for c. m. of the 
clump. The axis x connects the c.m.'s of two clumps 
when the distance between them is the shortest. One 
approximately has 

v x - v x ~ ^-AZ = ^-rcosip, (72) 
ol ol 

where ip is the polar angle in spherical coordinates. 

For nearly straightforward propagation, the angle be- 
tween Vj-d and the line connecting c.m.'s of the clumps 
evolves as 



dt 



Vrcl 2 
—j—COS i 



(73) 
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where u ro i — \2Vh- Changing variable t to (f> in the 
Newton equation one gets 



Vrell 



COS ( 



and after integration of this equation: 



2GM' 



V IC \R 



j9(y), 



(74) 



(75) 



where y = l/R', 



g(y)- 



-\ y>i; 
-(i-y 2 ) 1/2 ]/y, y<h 



and 2F1 is the hypergeometric function. 

The rate of internal energy growth due to collisions 
with all other clumps is 



E= / 2nlv TeX dl I dM'^(M',t)AE. 



(77) 



After simple calculations, one obtains 



E = 



4tt(3 - /3) G MR? 

3(5 - /3) W ro l 



dM'M' 2 ^{M',t) 



A I 



' A 






.i?' 2 " 




Vi?' 2 



where 



dy y 



0.11 



(78) 



(79) 



for x' c = 0.1 and dependence of A on a;' c is very weak. 

As the second process of tidal destruction we shall con- 
sider the interaction with the common gravitational po- 
tential of a host clump. Energy gain per mass unit at a 
distance r from c. m. of the considered small clump m 
during one periastron passage [5(j is 



(E P ) = 



GM h 
~R^'' 



Rh 

R P 



Xccc(e) A(ujt), 



(80) 



where e is the eccentricity, the function Xccc presented 
in the adiabatic correction A(x) = (1 + x 2 )~ 7 , 

7 ~ 2.5 — 3 and R p is the periastron separation. The 
energy gain of a clump during one period T or b is AE = 
j (E p ) p int (r)d 3 r , and the rate of energy growth is 



E 



2AE 



(81) 



The rate of energy growth due to the both of aforemen- 
tioned processes is given by the sum of (f75|l and (|HT|l . 



Using the distribution (|69|l in the integral of l|78(l and 
the total energy of a clump given by Eq. (|70|l we find: 



' - - 2t7 1 /^ /2 ^ 3 m("+ 3 )/ 2 mr 3( " +3)/4 , 



T{m,m h ,v,v h ) ~ E ~ A ^ V 
where 



(82) 



2VV/2(5-2/3) 



+ 



3(5-/3) 



n + 3 



) .lUT)\,. r , (< ) 



4^1/2 



n + 1 



2A 



n + 5 
n + 9 



(83) 



The first term in the square brackets describes the inter- 
action with common gravitational field of the host clump, 
while the second term — "collisions" with small clumps 
inside the host clump. Usually the former is larger than 
the latter. For calculations we shall consider an aver- 
age orbit with R h /R p ~ 2 and T orb ~ 2R h /Vh, and put 
j4(wr)x ccc (e) ~ 1. If to neglect the tidal interactions 
with the small clumps (the second term in the square 
brackets in Eq. I|83(l ') and to use j3 — 1.8, one has /i ~ 9.6. 
The dependence of our final result (mass function of the 
clumps) on fi is weak, approximately as fi^ 1 ^ 3 . 



C. Survival probability in the hierarchical model 

A small clump with mass m during its life time can be 
a constituent part of many host clumps of successively 
larger masses ml and virial velocities V . After tidal dis- 
ruption of the lightest host, a small clump automatically 
becomes a constituent part of the heavier host etc. Tran- 
sition of a small clump from one host to another occurs 
almost continuously in time up to formation of a host 
where tidal destruction becomes inefficient. A fraction of 
small clumps with mass m escaping the tidal destruction 
is given by e~ J , with 



At 



•f T(m, to', v, i>') ' 



(84) 



where summation goes over the intermediate big clumps 
which successively host the small clump to, and At ap- 
proximately equals to the difference of formation times 
of two successive hosts. 

Let us introduce the notation: mi is the mass of the 
first (lightest) host clump which contains the considered 
light clump m, and m n is the last such object, e. g. the 
Galactic halo. A formation epoch for the host clump to' 
is 



t f (m ,v)=t A l- 



l + z' 



1 \ -3/2 



+ Z\ 



tATO^/V- 3 / 2 . 



(85) 

Note that J does not depend on t\ since it enters linearly 
in tf and in T as seen from Eq. (|82|) . 
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The first host gives a major contribution to the clump 
destruction, especially if its mass mi is close to m and 
v\ ~ v. For the considered clump m and for the first two 
hosts cr eq (M) ~ cr cq (Mi) ~ er q(^2) because of a eq (M) 
depends weakly on M according to l|64l) . The considered 
clump m and the first two hosts differ mainly by masses 
and by values of v (y > v\ > v^). It justifies the following 
simplification: we consider all hosts beginning from the 
second one as the typical objects (i/j = 1 for i = 2, 3, . . .) 
and separating the first term in the sum (|84p. substitute 
the remaining sum by the integral. 



J 



tf{m-2,v 2 = 1) -t/(mi,i/i) 



+ 



T(m, mi, i/, z/i) 
1 



T(m, m', v, v' = 1) cim' 



df / (m',i/ = l) dm# _ (86) 



FIG. 2: Function from Eq. 1901 obtained by the numeri- 
cal integration of l|88jl . This curve is valid with good accuracy 
for all j3 from the interval 1 < f3 < 2. 



We may change now the lower limit m% to mi in the 
above integral because it depends on mi weakly, only 
through m^ 3 " 2 with (n+3)/2 <C 1, where n is given by 
Eq. i|64|) . I n addition we may put the upper limit m n — > 
oo without loss of accuracy. Inserting Eq. (JH21 into 
we obtain finally the following approximate expression 
for J: 

(n+3)/2 

J 



//.;/ 



9/2 


/ m 




\mi 
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m \ ' 



(n+3)/2 



(87) 



where the step-function Q(x) = 1 for x > and = 
for a; < 0. In the case of v\ = 1 the formation moment 
of the first host almost coincides with the formation mo- 
ment of the second one, for which vi = 1 in the used 
approach. 

The differential fraction of mass in the form of clumps 
which escape the tidal destruction in the hierarchical ob- 
jects can be found as 

i(n,v)—dv = dmdis(2/ir)e- u2/2 f dv x e~ v *l 2 
m J 



I dti 



d 2 F(m,t 1 ) 



dm dt\ 



e- J , (88) 



where the variable mi in J is connected with t\ and v\ 
according to Eq. Q85JI. F(m,t) is the mass fraction of 
unconfined clumps with masses smaller than m at the 
moment t, which according to |40| is given by 



F(m, t) = erf ■ 



6 C 



V2a cq (M)D g (t) 



(89) 



Here erf (x) is the error- integral and D g (t) is the growth- 
factor. After numerical calculations for (3 = 1.8 we finally 
obtain: 



e(n,^)^(27r)- 1 /2 e - 2 /2 (n + 3)yM5 



where the function y(v) is plotted in the Fig.0 It weakly 
depends on (3. 

This is the main result of this section. We shall refer 
to £ as to the clump survival probability. 

Our final result, the distribution of surviving clumps 
over their masses M and v, which for the fixed M deter- 
mines the clump radius R, is given by £(n, v) dudM/M, 
where £ depends very weakly on M only through the 
weak dependence of n(M) (see Eq. (16411 ). For most nu- 
merous clumps with v = 1 and for the unit intervals dv 
and dhxM, £ has a meaning of fraction of DM mass in 
the form of clumps relative to free DM particles, as it was 
introduced in Eq. (|10|l . By integrating over we obtain 



£ int ~0.01(n + 3). 



(91) 



(90) 



This means that, for different n, about 0.1 — 0.5% of 
clumps survive the stage of tidal destruction in each log- 
arithmic mass interval A In M ~ 1. 
Several remarks are in order. 

The physical meaning of the survived clump distribu- 
tion £(n, v)dvdM /M is different from that for unconfined 
(free) clumps given by the Press-Schechter mass func- 
tion dF/dm. Surviving clump distribution implies that 
each DM particle belongs simultaneously to several host 
clumps put into each other, and by this reason the inte- 
gral J £dm/m is divergent. 

To calculate £ we need to know the power-law index n 
in the perturbation spectrum, which can be taken as n e g 
from Eq. I|fci4fl for given <j eq (M). To find the latter, the 
primeval (e. g. inflation) power spectrum of fluctuation is 
needed. The simplest inflation models give P p (k) oc k np 
with n p m 1. The analysis of the WMAP measurement of 
CMB anisotropy |5l| gives power law spectrum with n p — 
0.99 ± 0.04 in good agreement with n p — 1. However, 
when data from 2dF galaxy power spectrum and Ly-a 
are included in analysis, the best-fit favors a mild tilt, 
n p = 0.96 ±0.02. 

The variance a eq (M) = U( ) from Eq. H44fl in the small 
scale range is found in [5^] (see also |53^ . We present 



13 



this result as 



where 



o-eq(M) 



2 x 1(T 4 



hi 



icq 



3/2 / fc \ ("p-1)/2 



8 ^ 2x1 q3.7(™ p -1)-3 



/ M 

1 -°-° 61g (^ 



(l-n p )/6 



3/2 



(92) 



where fc oq and k^o correspond to mass inside the cosmo- 
logical horizon at the moments t eq and to, respectively, 
and 

/ S (0 A ) = 1.04 - 0.82f2 A + 2f2 A (93) 

according to [H^- We used above the relation k oc M^ 1 / 3 
and the values 



M eq = 1.5x10 



i!) 



n m 'h^g, M ho = 6xlO bb n m h- i g. (94) 



Using the power-law spectrum of fluctuations down 
to the small scales, while normalization by the CMB 
anisotropy is performed at large scale, implies an extrap- 
olation of the spectrum by many orders of magnitudes. 
Such extrapolation is justified only by the confidence to 
inflation models which predict the power-law spectrum 
valid for many orders of mass magnitudes. 

It is interesting to note that the differential num- 
ber density of clumps in the Galactic halo n{M) dM oc 
dM/M 2 , obtained from Eq. I|9U|I . is very close (includ- 
ing the normalization coefficient) to that obtained in 
the numerical simulations for clumps with large masses 
M > 10 8 M Q [n(M)dM oc dM/M 1 - 9 [|]. Strictly speak- 
ing our calculations are not valid for clumps with these 
masses, because of their destruction in the halo up to 
present epoch to and accretion of new clumps to the halo. 
Nevertheless, for the small interval of masses where the 
power-law spectrum can be used as a rather good ap- 
proximation, our approach can be roughly valid. 



VI. NUMERICAL RESULTS 

In this section we shall present the numerical results 
of our calculations: enhancement factor for annihilation 
signal, the distribution of DM clumps over their masses 
M and radii R, and the distribution of clumps in the 
galactic halo. 

Using Eqs. © and (JZJ) we find the enhancement of 
the annihilation signal due to dumpiness of the halo as 
generalization of Eq. I|1(J|) for the clumps distributed over 
M and R: 



V{Mnun,n p )=l + - 



1 



dp 



Pdm 



'dM 
W 



S\j3,x c (v)]£p int (M, v), 
(95) 



£ = £(n,i/) 



(96) 



is defined by Eq. I|90(l . an effective spectrum index 
n(M,n p ) is calculated from Eq. (|64|l for primeval (infla- 
tion) spectrum index n p , with a c s (M) taken in the form 
(EHJ); pint is given by Eq. (IB1)|) and v min ~ 0.55. Function 
S is taken in the form Q , which corresponds to the den- 
sity profile and we used Eq. I|60|) for a clump virial 
radius R. 

The major part of the survived clumps are formed from 
fluctuations with a mean value of the peak height 



dv ve v l 2 y{v) 
dve-^ l2 y(v) 



1.6. 



(97) 



Meanwhile the main contribution to the enhancement of 
the annihilation signal (|95|l comes from the clumps with 
an effective value of the peak height: 

/ \ dvvS£p- mt dM/M 

(^)ann ^ 4^7 ^ 2.5 (98) 

dv S£p int dM/M 

for (3 — 1.8. The clumps with v ~ 2.5 have x c ~ 0.05. 
For the Galactic halo we use the NFW density profile 



Pdm(1) 



IX) 



(l/L) (1 + l/L)) 



2 ! 



(99) 



with L — 45 kpc according to |&4|, and po fixed by the 
local density value pdm^o) = 0.3 GeV cm" 3 . With 
these parameters the halo mass within the virial radius 
of 100 kpc is 10 12 Af o . Eq. © gives p H = 1.02p H (r Q ), 
i. e. these values practically coincide. 

The values of global enhancement r/(M m i n , n p ) as given 
by Eq. I|95ll are displayed in Figs. for different values 
of M m i n , (3 and n p . As a representative example consider 
the clump with the Gurevich-Zybin density profile 
with (3 = 1.8 (see Fig. [3J : numerically i] = 5 for M min = 
2x 1O _8 M0 and n p = 1.0. It strongly increases at smaller 
M m in and larger n p . For example, for n p = 1.1 and 
n p = 1.2 at the same M m i n = 2 x 10" 8 M Q , enhancement 
becomes tremendously large, rj = 130 and rj = 4 x 10 3 , 
respectively. 

Our approach is based on the hierarchical clustering 
model in which smaller mass objects are formed earlier 
than the larger ones, i. e. a eq (M) diminishes with the 
growing of M. This condition is satisfied for objects 
with mass M > M min ~ 2 x 1O~ 8 M only if the pri- 
mordial power spectrum has the value of the power in- 
dex n p > 0.84. As seen from Figs. |3HHI m this case the 
enhancement of the annihilation signal in fact is absent, 
i] ~ 1, for n p < 0.9. 
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FIG. 3: The global enhancement rj of the annihilation signal 
from Eq. I|950 as a function of the minimal clump mass M m i n , 
for clump density profile with index = 1 and for different 
indices n p of primeval perturbation spectrum. The curves are 
marked by the values of n p . 
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FIG. 4: The same as Fig.glbut for j3 = 1.5. 



We discussed above the global enhancement of the 
annihilation signal. In fact the enhancement varies in 
the different directions relative to the Galactic Center 
(GC). It can be easily seen from Eqs. © and Q, which 
show that while / c i is proportional to n c \(l) (and thus to 
/3dm (0 everywhere except the core), 4om is proportional 
to Pdm(0- It implies that relative contribution of /hom 
increases in the directions close to GC where Pdm(0 is 
larger. This effect is further enhanced by the destruc- 
tion of the clumps in the core around the GC. Our nu- 
merical calculations confirm this expectation: the ratio 
of enhancement in the directions to GC and Anticenter 
is 0.2 for NFW density profile and for the core radius 
3 kpc. Therefore, the presence of the clumps diminishes 
the anisotropy of the annihilation signal, caused by the 
density profile of DM in the halo. 

Our results suffer from uncertainties in input param- 
eters. As was mentioned above /5dm — Pdm(^©)- It re- 
mains approximately true not only for the NFW den- 
sity profile but also for other profiles discussed in the 
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FIG. 5: The same as Fig.^Jbut for /3 = 1.8. 

literature. For example, for isothermal profile with the 
core radius 10 kpc, /5dm = 0.65pDM(f©)- Another 
uncertainty in the value of rj is imposed by value of 
Pdm(^o)- According to different estimates pdm^o) = 
(0.2 — 0.6) GeV cm" 3 . The corresponding uncertainty in 
rj is given by factor 0.5 — 1.5. 

For illustration we shall describe numerically the prop- 
erties of clumps which give the main contribution to an- 
nihilation rate. Basically, they are those with M ~ M m [ n 
and v ~ (v). The r.m.s. fluctuation values a eq (M) for 
clumps with minimal mass Af m ; n ~ 2 x 10 _8 Mq and for 
rip = 1 and 1.2 are er cq = 0.015 and 0.14, respectively, ac- 
cording to Eq. ■ From Eq. the effective value of 
v = 5 cq /a cq is (v) ~ 2.5. From Eqs. ® , @3 it follows 
that clumps with this v are characterized by the density 
and radius p int ~ 2 x 10~ 22 g cm -3 , R ~ 3.6 x 10 15 cm 
and pi n t ~ 2x 10~ 19 g cm -3 , R ~ 3.7 x 10 14 cm for n p — 1 
and 1.2, respectively. The part of Galactic halo mass in 
the form of these clumps is of the order of £i n t ~ 0.002 
according to Eq. 1)91(1 . A mean number density of the 
clumps in the halo is ~ 25 pc -3 . 

We have given above the characteristic values for 
clumps with dominant contribution to the annihilation 
signal. The general distribution of clumps in the Galactic 
halo can be readily calculated numerically from Eq. 19011 , 
changing (for given M) the distribution over v by that 
over dR: 

n c i(M,R)dlnMd\nR = PDM n ^ Q - £(M, v) d In M dv. 

( 10 °) 

Note that definition of the clump number density n c \ here 
does not coincides with the similar one in Eq. JSJ where 
dN is given per dR and dM. 

The distribution Mn c i(M, R) is presented in Fig. [fj] as 
function of R for different M and for distance 8.5 kpc 
from the Galactic center. 

The radius of a clump R in most general case is deter- 
mined by M and v. Due to weak dependence of a eq (M) 
on M, the radius of the clump R(M) with sufficiently 
good accuracy is proportional to M 1 / 3 . From Eqs. I|66|) 
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FIG. 6: The mass density of clumps in the Galactic halo 
M n c i(M,R) in units M / kpc 3 , from Eq. I100L as function 
of their radius R at the distance 8.5 kpc from the Galactic 
center for n p — 1.0. The curves are labeled by values of clump 
masses in Mq 



and I|b7fl we have 



R~ 1.5 x 10 16 



M 



10- 6 M ( 



1/3 



v 

2J3 



cm, 



(101) 



for n p — 1.0. For n p — 1.1 the numerical factor in the 
Eq. ifTUTl) is 5.5 x 10 15 cm. 

How clumps are distributed in the Galactic halo? One 
may expect that this distribution is the same as distribu- 
tion of the free DM particles in the halo. This is true for 
large distances I from Galactic Center, while at small I 
the tidal interaction with stars results in the destruction 
of small clumps and in the formation of the core with 
radius L c in the clump distribution. The destruction 
of a clump propagating in the space filled by stars has 
been studied in |9(. The destruction is important only 
for clumps inside the bulge, i.e. at the distance I < 3 kpc 
from the Galactic Center. Clumps outside the bulge at 
the distances 3 kpc < I < 10 kpc can interact tidally 
with stars in the disk. But time of crossing the disk is 
very small (in comparison with orbital period) and this 
process is not important. 

The number density of clumps outside the bulge is pro- 
portional to the halo density, e. g. to £/9dm(0 in the case 
of the NFW distribution given by Eq. (|99l) . or can be 
obtained from (|99|l and Fig. by simple scaling. 

The density distribution for stars in the bulge accord- 
ing to (5f| is given by 



Ps(l) 



p(l/f) 



-1.8 



I < 



p(l/f)- 3 , I > r, 



(102) 



where p — 1.8Af Q /pc 3 and f = 800 pc. From Eq. (46) 
of by substituting Eq. we obtain that inside 

the bulge (/ < 3 kpc), the clumps with M < 10~ 4 M o 
are destructed during the Hubble time. Thus, for these 
masses the core radius L c coincides with the size of the 
bulge Lbuige ^ 3 kpc. 



log(M/M Q ) 

FIG. 7: The radius L c of the Galactic core (in kpc) in the 
distribution of clumps with masses M for n p = 1.0 and n v = 
1.1 in the case v = 2.5. 



Clumps with M > 10- 4 M Q are destructed during the 
Hubble time within distances from the Galactic Center 
shown in Fig.0for n p = 1.0, n p = 1.1 and v = 2.5. This 
distance defines the radius of the core L c for clumps of 
the given mass M. 

Our calculations for enhancement of annihilation sig- 
nal disagree with those in 0, |27], |28j . 

In 0| the singularity in the Galactic Center is cut 
at very small core radius, which results in too strong 
annihilation signal. According to our calculations the 
radius of the core is much larger, and the distribution of 
the clumps in the halo also has a core. 

In [2^.l28| the large enhancement of the signal is found 
for heavy clumps with M > 10 6 M Q [13 and M > 1O 2 M 
[2c| . If it were true, the total signal from clumps with 
AI > Af m i n would be too large. Too small core radius 
was used in these calculations, too. 



VII. CONCLUSIONS 

We have calculated the number density of the small- 
scale clumps in the Galactic halo and their distribution 
over masses M, radii R and distances to the Galactic 
Center in the framework of the standard cosmological 
model with the primeval density perturbation P(k) oc 
k np taken from the inflation models with n p ~ 1 (the 
Harrison-Zeldovich spectrum). The most important el- 
ement of our calculations is inclusion of the tidal inter- 
actions, which result in the formation of the clump core 
and destruction of small-scale clumps. 

We consider most conservative case of the Gaussian 
adiabatic fluctuations which enter the non-linear stage of 
evolution, t n \, at the matter-dominated epoch t n \ > t cq , 
where t oq is the moment of equality. The time of small- 
scale clump formation tf for a clump with mass M is 
given by two equation: formation criterion 8{M, tf) = S c 
and by height of peak density of a fluctuation in units of 
dispersion v = S eq /a eq (M), taken at the epoch t eq [see 
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Eqs. i|48|) and (|62() for explanation and notation]. All 
processes we are interested in, take place at t > t cq at the 
stage of non-linear evolution. The growth of fluctuations 
in the non-linear regime we study in the framework of 
the Press-Schechter theory of hierarchical clustering with 
the tidal interactions included as the new element. The 
picture of hierarchical clustering and clump destruction 
can be described in the following way. The clumps of 
the minimal mass are formed first. A clump of larger 
mass, which hosts the smaller clumps, is formed later. A 
bigger clump, which includes the considered hosts with 
their content, is formed further later etc. The clumps are 
destroyed in tidal interactions with other small clumps 
and by gravitational field of a host clump, with the former 
process being subdominant. The calculated mass density 
of the survival clumps, £(M, v)dvdM /M , with mass M 
and v, is given by Eq. I|90[) . with the survival probability 
being typically £ ~ 0.001 — 0.005. The clump number 
density in the Galactic halo n c \(M, R) for the clumps 
with mass M and radius R is shown in Fig. [S] These 
clumps are distributed in the Galactic halo as function of 
a distance to the Galactic Center I. At large distance the 
distribution must be the same as found in the numerical 
simulations (e. g. the NFW profile). At small distance 
there is a core produced by tidal interaction of the clumps 
with the stars in the bulge. The radius of the core, L c is 
given in Fig. at M > 10 _4 Mq and it is equal to the 
radius of the bulge L c ~ 3 kpc for smaller clump masses. 

The mass spectrum of the clumps is characterized by 
a cutoff at M m j n . Its value depends on the properties 
of the DM particle, and thus it is model dependent. 
The existing calculations of M m - m differ drastically: from 
M min ~ 1O- 12 M to M min ~ 1O- 7 M 0. 

Cold Dark Matter particles at high temperature T > 
Tf ~ 0.05m x are in thermodynamical (chemical) equilib- 
rium with cosmic plasma, when their number density is 
determined by temperature. After decoupling at t > tt 
and T < Tf, the DM particles remain for some time in 
the kinetic equilibrium with plasma, when temperature 
of CDM particles T x is equal to temperature of plasma 
T, but number density n x is not Planckian any more. At 
this stage the CDM particles are not perfectly coupled to 
the cosmic plasma. Collisions between a CDM particle 
and fast particles of ambient plasma result in exchange 
of momenta and a CDM particle diffuses in a space. Due 
to diffusion, the DM particles leak from the small-scale 
fluctuations, and thus their distribution obtain a cutoff 
at minimal mass Mo- The diffusion coefficient is deter- 
mined by elastic scattering of DM particles off the plasma 
particles. Our calculations made for neutralino, for which 
we have chosen the pure bino state, give 

M D = 4.3 x 1Q-13( m - ) - 15/S ( JL_) 3 M , 
VIOOGcV/ ^lTcVy °' 

(103) 

where m x is the neutralino mass and M is (approxi- 
mately) the mass of sneutrino and selectron, which are 
assumed to be equal. The functional dependence of 



Eq. H25f) and numerical value of Eq. (|26|l obtained in the 
diffusion approximation coincide with the corresponding 
results obtained by different method in |4^ . 

When the energy relaxation time for DM particles, r re i, 
becomes larger than the Hubble time i/ _1 (i), DM parti- 
cles get out of kinetic equilibrium. This conditions deter- 
mines the time of kinetic decoupling td- At t > td CDM 
matter particles are moving in free streaming regime and 
all fluctuations on the free-streaming scale A/ s and less 
are washed out. In contrast to 43], we have calculated 
the free-streaming length A f s taking into account the dis- 
tribution of neutralino (bino) velocities over absolute val- 
ues and angles from radial directions. The cosmological 
expansion in the vicinity of i oq is taken exactly, without 
usual step-function approximation. Our value of M m { n 
due to free-streaming effect is 

M ->°= i - 5 * i °- 8 (io^r /! (Tiv) 

(104) 

see Eq. I|38(l for more parameters involved. 

When normalized to the same masses of neutralino and 
sleptons, our value of M m ; n coincides only by order of 
magnitude with |43j. 

The evolution of a density fluctuation in the non-linear 
regime results in the density profile of a clump. The 
analytic theory of this phenomenon was developed by 
Gurevich and Zybin (for a review see |9j ) , for the numer- 
ical simulations see The initial single-stream flow 
leads to formation of initial singularity. In contrast to 
energy-dissipating matter (e. g., baryons), in the flow of 
non-dissipative matter the multistream instability devel- 
ops [9j, when at one point several streams with different 
radial velocities exist. The surfaces with different num- 
ber of streams are separated by caustics, which number 
increases rapidly towards the center. The matter is grav- 
itationally captured in a such structure. The density sin- 
gularity is produced in the center, unless the additional 
phenomena arc included in the consideration. As such 
the interaction with the damped mode |ij and annihila- 
tion of DM particles [2j| were previously studied. Wc 
have demonstrated here that tidal forces due to external 
gravitational field cause the deflection of DM particles 
from radial motion, and prevent thus formation of singu- 
larity. The produced core has a radius R c given in the 
approximate form as 

x c = % ~ 0.3^ 2 , (105) 
R 

(see Eq. (|60H for the exact expression and the discussion 
afterwards). This radius is much bigger than those ob- 
tained in [alii- 

The majority of clumps are formed from v ~ 1 peaks, 
while the survived clumps are characterized on average 
by v ~ 1.6. The clumps which give the dominant contri- 
bution to the annihilation signal have v ~ 2.5. 

In spite of small surviving probability £^0.1 — 0.5%, 
clumps in most cases provide the dominant contribution 
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to the annihilation rate in the halo. The enhancement 
of the annihilation signal can be characterized by ratio 
V = (-^ci + ^hom)/^hom, where I c i is the annihilation signal 
from the clumps, and /hom - from homogeneously dis- 
tributed DM particles with the NFW density profile in 
the Galactic halo. The main contribution to r\ is given 
by v ~ 2.5 and M ~ M m i n . The signal enhancement rj is 
shown numerically in Figs. [3H3 One can see that practi- 
cally for all allowed values of primeval perturbation spec- 
trum index n p > 1.0 the annihilation signal from clumps 
gives the dominant contribution. This result does not 
depend on the properties of DM particles. 

The observations favor the spectrum index n p = 1.0 
pHj • The enhancement of the annihilation signal for this 
value of n p is described by factor 2 - 5 for different (3 with 
uncertainties due to values of M m i n and other parame- 
ters. 

The clumps which give the dominant contribution to 
the annihilation signal have approximately the following 
properties in the case n p = 1: The mass M ~ M m j n and 
v ~ 2.5, the radius R ~ 3.6 x 10 15 cm and the radius of 
the core R c ~ 1.8 x 10 14 cm, the mean internal density 
of the clump p lnt ~ 2.5 x 10~ 22 g cm , the fraction of 
the halo mass in the form of these clumps £i„t ~ 0.002, 
and the mean number density of these clumps in the halo 
n c i ~ 25 pc~ 3 . 

Recently HEAT collaboration detected excessive flux 
of cosmic ray positrons at energy E ~ 10 GeV j5||. Ac- 
cording to |57|. if this positron flux is produced by an- 
nihilation of neutralinos the enhancement factor of order 
of 30 is needed. The calculations presented here show 
that such enhancement can be reached in the considered 
scenario in case of extreme combination of parameters. 
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APPENDIX A: CROSS SECTIONS OF 
NEUTRALINO SCATTERING OFF ELECTRONS 
AND NEUTRINOS 

As neutralino we shall consider here a pure bino (x — 
B). The Lagrangian for interaction of bino with left and 
right components of a fermion / can be written (see e. g. 
[Hi HI) as 

C ffx = -gV2tzn6 w (e f -T 3 L f )fP RX fL 

+gV2 t&n0we f fPLXfR, (Al) 

where g is SU(2) coupling constant, 9 W is the Wein- 
berg angle (sin 2 #w = 0.231), et is electric charge of the 



fermion / in the units of electron charge, is the pro- 
jection of weak isospin for f L , P R = 1/2(1 + 75) is a 
projection operator which cuts the left component from 
the operator / in Eq. (|A1|) ; Jl is the left sfermion. The 



first term in the Lagrangian i|Al|l is C 



C 



the second 



~ , » . When f is ultra-relativistic in the frame where 

JrJrX. j 

neutralino is at rest, there is no interference for scatter- 
ing of the left and right components of the fermions (the 
interference terms are proportional to rrif). Therefore, 
we shall calculate fx cross-section for left Jl and right 
fn fermions separately. 

Scattering of the left fermion with e/ = —1 and T3/ = 
— 1/2 (e. g. e, //, t) off bino are described by the two 
diagrams in s- and it-channels as shown below 




f L k 2 k x 



X\ Vi Pi ■ 



h 




P2 



1 

j , 



x h 
u 



The standard calculations for matrix elements give for 
|M| 2 = |M S | 2 + \M U \ 2 + 2Re(M s M„*): 



1 



■2 _ ' ,,( A-|/;i )\k,p,) 

(S-TO 2 ) 2 



\M s y = -(fltanflw) 



(pife 2 )(fclP2) 



m 2 L ) 2 



mm:, 



-(gtan^ 



rrPikik?) 



(s-m 2 )(u-m 2 )' 



(A2) 



(A3) 



(A4) 



Cross-section for the Jl + X ^ Jl + X scattering at angle 
012 in the system where neutralino is at rest is given by 



da c \ 
~dtt 



1 



64tt 2 s 



M\ 



w 2 (l + cos6»i 2 ) 



8 cos 4 9 



w 



(m 2 



rh 2 L ) 2 



(A5) 

where u> ^> mf is energy of Jl in the system where neu- 
tralino is at rest, m x is the neutralino mass and is 
the mass of the left sfermion. 

Let us consider now fn + X — * fn + X scattering de- 
scribed by the second term in r.h.s. of Eq. IjAljl . The dia- 
grams are identical to that in the figure after substituting 
/l — * /r and /r, — * /r. Since traces do not change when 
Pl Pr, the expressions (|A2|) - (|A4|) remain the same, 
changing only due to coupling constant which increases 
twice (see Eq. l|Al|l ). Therefore, we obtain 



Irx V fl "//LX 

In this paper we are interested in v + x ~ ¥ v + X an d 
e + x ~ > e + X scattering. In the former case the cross- 
section is given by Eq. I|A5(1 , and in the latter case - by the 
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sum of f L + x -> h + X and fx. + X -> fit + X scattering, 
i.e. it is by factor 17 larger than the cross-section l|A5[) . 



APPENDIX B: KINETIC EQUATION 

In this Appendix we shall study the stage of kinetic 
equilibrium and the stage after its breaking in the com- 
mon formalism of kinetic equation similar to |42^ and 
using approach of |4J| . We shall confirm in this way the 
results of Sec. lIIII and clarify the difference in calculations 
of M min . 

Following [44j we introduce the neutralino distribution 
function f(x,p,t) over comoving coordinates x and mo- 
menta p — ma?x (with this definition momentum is con- 
stant for freely moving particles). The neutralino density 
is 



p(x,t) = ^ J d 3 pf(x,p,t)=p x {t)(l + 5(x,t)). 



(Bl) 



The kinetic equation with the collision term of the 
Fokker-Planck type 0] can be written as 



df | Pl df m 8<j> df 

dt ma? dxi dxi dpi 



p{) d Pi \mTa? J + d Pi 



(B2) 

where <fi is the gravitational potential, which can be ne- 
glected at the considered epoch t < t cq , T(t) is the tem- 
perature of the ambient plasma given by Eq. O, and 
D p (t) is the diffusion coefficient in the momentum space. 
According to 0] 

D P (t) = ^-JmJ dunoitu) (jgf) (Sp) 2 . (B3) 

The number 40 in Eq. I|B3|I comes from the counting of 
degrees of freedom in neutralino-fcrmion scattering as in 
the Sec. EH 

The equation l|B2|) with the diffusion coefficient l|B3|) 
coincides with Eq. (16) from (4^ except the numerical 
factor in D p which is of order of unity. 



1. Kinetic decoupling 

Let us consider an exit of neutralinos from the kinetic 
equilibrium (decoupling) in the homogeneous universe, 
when d/dxi terms in Eq. (|B2() can be neglected. The 
temperature of neutralino gas T x is defined as 



PiPjfd 3 p = p x a 5 T x (t). 



(B4) 



Multiplying Eq. I|B2|) by PiPj and integrating it over d 3 p 
one obtains 



dT x 
dt 



2 d T 2D p (t) 
a x ma? 



1 - 



T x (t) 
T(t) 



= 0, 



(B5) 



The initial condition for Eq. I|B5I) can be chosen at the 
moment of freezing t — t j as in |42| , or more conveniently 
at any t% from interval i/ < ti <C td, as T x (ti) = T(ti), 
where T is the temperature of ambient plasma. Solu- 
tion of Eq. i|B5(l (see below) gives transition of ratio 
r(i) = T x (t)/T(t) from r = 1 to < 1 within some time 
interval, determined by r^. Any value of t in this interval 
can be taken as definition of decoupling time td- Eq. JB5I) 
and its solution can be simplified using the dimension- 
less time t = t/td- Characteristic time td is naturally 
emerged from dimension parameters entering the diffu- 
sion coefficient, and up to numerical coefficient it coin- 
cides with td determined in Sec. IIHI The transition time 
interval fixes this numerical coefficient with some uncer- 
tainty, and we obtain indeed td (and hence Td = T(td)) 
approximately equal to those given by Eqs. i|19fl and l|2(J|) 
in Sec. I11II The solution of Eq. i|B5() in terms of r = t/td 
is given by 



V/ 4r2 d^x^e 1 ^ 2 




(B6) 

The asymptotic forms of solution (|B6J| are given by 
T x /T d = t- 1 ' 2 for r < 1 and T x /T d = r- 1 r(3/4)/2 1 /2 
for r ^> 1 as it must be. From solution (|B6|) it is seen that 
transition from kinetic equilibrium of neutralino with rel- 
ativistic fermions to the non-equilibrium regime proceeds 
very fast. By this reason our consideration of diffusion 
and free streaming independently in Sec. IIIII is well jus- 
tified. 



2. Diffusion 

Consider Eq. l)B2jl before kinetic decoupling, t ; <C td- 
One can find the first two moments by integrating first 
time Eq. I|B2|) over d 3 p and second time over pid 3 p. In- 
serting the first of the obtained equation into the second 
one we obtain the following equation for the Fourier com- 
ponents: 



<^5 

d 2 t 



a 85 



D p {t)- 



1 85 



Px 



' adt ' py "'mTa? dt 
The r.h.s. of the IB 711 has a tensor form 



-f— / p l Pjfd i p. 
a'm J 

(B7) 



f 



p x a 7 m 



PiPjfd 3 p = E5 tj + Fkikj 



(B8) 



where the isotropic part E — T x 5k/a 2 m for any r, while 
F depends on time t. In the limit t « 1 we may put 
F = and neglect the first and second terms in (|B7|I . 
The resultant equation coincides with diffusion equation 
(|24Jl with the same diffusion coefficient Ij23(l and has the 
same solution. 

In ^3] only this diffusion limit of the general kinetic 
equation l|B2|l has been considered. 



19 



3. Free streaming 

In the limiting case t 1, i.e. after decoupling, the 
B2I has a simple form 



df | Pi df 
dt ma 2 dx,; 



= 0, 



(B9) 



according to (|B6(1 . kinetic decoupling proceeds very fast. 
Integrating l|B10[) over d 3 p with initial condition 



f{t d ) = {2nT d ma 2 d )- 3 ' 2 e^ 



2T d ma d 



(Bll) 



with the solution 



/ oc exp 



ma d 



g(t) 



(BIO) 



where g(t) is the same function as J2HJ- The solution 
(IBlOf) is valid with a good accuracy also at r > 1, because 



one obtains 



and then Eqs. |SJ|, JBU and J3SJ| from Sec.|ITT| 



(B12) 
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